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I. INTRODUCTION 



In writing this article, I have two main objectives. The first is to chart the evolution of 
the analytic approach to Quantum Chromodynamics (QCD) perturbation theory, focusing 
on the key turning points in its theoretical development. The second is to demonstrate how 
this approach works in terms of selected applications. Of course, I am forced to leave out 
many interesting applications and aspects of the approach, which in turn means that my 
presentation will not be complete. I would refer the reader to the original references for 
further reading and in order to study the subject in more technical detail. The sorts of ideas 
and issues to be discussed in this presentation include the following benchmarks: 

• "Ultraviolet (UV) freedom" , "infrared (IR) slavery" , and Landau singularity. 

• Renormalizability, analyticity (in Q 2 and g 2 ), causality, and summability. 

• First remedies in the infrared: Color saturation, effective gluon mass. 

• Screening of the Landau singularities via Sudakov form factors. 

• Shirkov-Solovtsov analytic coupling — Euclidean and Minkowski space. 

• From a recipe to a paradigm: Analytic Perturbation Theory (APT). 

• More scales, more riddles: Logarithms of the factorization (evolution) scale — Naive 
and Maximal analytization. 

• Generalization of the analyticity requirement: From the coupling to the whole hadronic 
amplitude and its spectral density. 

• Creation of Fractional APT (FAPT) in the spacelike and timelike regions. 

• Summation of perturbation series, Sudakov gluon resummation (exponentiation vs. 
analytization), power corrections. 

II. ULTRAVIOLET FREEDOM AND INFRARED CONFINEMENT 

Two key properties define QCD as the unbroken SU(N C = 3) (with N c denoting the 
number of colors) Yang-Mills theory describing strong interactions at the microscopic level: 
Asymptotic freedom and confinement. While the first property has been confirmed and was 
granted a Nobel prize in Physics in the year 2004, the other is still a matter of fact, with 
the underlying microscopic mechanism still at large. One claims that the reason for the 
non-observation of free quarks and gluons in nature is due to their color content. From this, 
one concludes that only color singlets are observable — hence, defacto color confinement. 

The main ingredient in justifying the UV freedom of QCD is that the renormalized QCD 
coupling becomes small at large momenta: g 2 — > as Q 2 — > oo. This implies that though 
initially the coupling and the momentum scale are independent parameters 1 they get linked 
to each other on account of the renormalizability of QCD. Hence, in the large-momentum 



Actually, there is no scale at all inside the QCD Lagrangian — provided one ignores current quark masses. 
A scale dependence is generated through dimensional transmutation. 
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region (or, equivalently, at short distances), QCD can be applied perturbatively, in analogy 
to QED in the sense of a weak-coupling expansion. A conventional way to define the running 
strong coupling in the MS scheme is given by (a s = g 2 /An) 
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which expresses it by means of the dimensional parameter A = Aqcd, whereas Nf represents 
the number of active flavors. For our considerations to follow it is more convenient to use 
the abbreviation L = ln(/z 2 /A 2 ) and recast the above expression in terms of the first two 
coefficients of the (3 function 
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which is an approximate solution of the Renormalization-Group (RG) equation 2 
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exhibiting the asymptotic-freedom property, i.e., a s ~ l/ln(/x 2 /A 2 ) — > as /i 2 /A 2 — > oo. 
The scale parameter A is not a generic parameter of the QCD Lagrangian; its definition 
is arbitrary. However, once defined, it becomes the fundamental constant of QCD to be 
determined by experiment. 

On the other hand, at fi 2 = A 2 the running strong coupling has a Landau singularity 
that spoils analyticity. At the one-loop level this singularity is a simple pole, whereas at 
higher loops the corresponding Landau singularity has a more complicated structure (e.g., 
a square- root singularity at two loops, et cetera). To restore analyticity in the infrared 
region < Q 2 < A 2 and ensure causality in the whole Q 2 plane, one has to remove the 
Landau singularity at any loop order of the perturbative expansion. This may become a 
serious problem, both theoretically and phenomenologically, for several hadronic processes — 
in particular for exclusive processes like hadron form factors. For example, in the case of the 
spacelike pion's electromagnetic form factor — a typical exclusive process — most experimental 
data are available in a momentum region, where repercussions of the Landau singularity can 



2 The first two coefficients f3o and j3\ are scheme independent. The next coefficients 02 and have also 
been analytically computed and confirmed by independent calculations at the level of 3 and 4 loops (see 
second paper in |l[ for references). 



3 



influence (actually contaminate) the perturbative calculation, as first shown and detailed in 
[2| and later further analyzed in Another way to describe this behavior is to use the 
language of information theory and consider a perturbatively calculable quantity, like the 
scaled pion factor vs. the momentum transfer Q 2 , as being a dynamic range of quality of 
the perturbative expansion. In this context, the accuracy of the QCD perturbative result 
(the signal) depends on the magnitude of the running coupling. The larger Q 2 , the smaller 
the coupling and, hence, the higher the precision of the (perturbative) calculation is. As 
one approaches the Landau singularity, the coupling increases rendering the perturbative 
expansion less and less precise, i.e., creating a kind of noise. When one reaches the Landau 
singularity, the "noise" (the uncertainty) becomes overwhelming and one is unable to extract 
any information (a clear "signal" ) from perturbation theory. The next few sections will show 
in more detail how this problem can be averted. 

Of course, one may interpret the breakdown of perturbation theory at low momenta 
~ A, as a "natural" consequence of the (unknown) confinement forces that prevail in this 
momentum region. The problem is that this breakdown is not universal but depends on the 
renormalization scheme adopted and the renormalization scale-fixing procedure used. For 
instance, adopting the Brodsky-Lepage-Mackenzie (BLM) commensurate-scale procedure, 
when including higher-order QCD perturbative corrections (see [2], 0] for further details), 
the BLM scale may become very low (compared to the original momentum scale) and deeply 
intruding into the "forbidden" nonperturbative region, where the running coupling blows 
up and perturbative QCD becomes invalid. An example: Using the BLM procedure in the 
calculation of the pion's electromagnetic form factor at NLO Q, the QCD running coupling 
may still "feel" the Landau singularity at momenta as large as 6 GeV 2 , a region close to 
the largest momenta probed by experiment at present for this reaction. Hence, the (scheme- 
dependent) optimization of the QCD perturbative expansion of typical hadronic quantities 
may depend in a crucial way on the IR behavior of the running coupling. It becomes 
clear that any attempt to apply an optimized renormalization prescription (scheme) in high 
fixed-order perturbative calculations should make sure that the result is not distorted by 
artificial logarithmic enhancements originating from the Landau singularity. Clearing away 
this obstruction is a real challenge. 

III. FIRST REMEDIES IN THE INFRARED 

The simplest way to avoid the Landau singularity in Euclidean space is to use a rigid IR 
cutoff on its value, say, a c s ntoS ~ 0.5 — 0.7, such as to render the perturbative expansion sound. 
A more sophisticated procedure is to assume that below some momentum scale /i IR > A color 
forces saturate, so that quarks and gluons are confined within color-singlet states. The IR 
cutoff scale /^ir, can be conceived of as an effective gluon mass m g generated by the dynamics 
of confinement, as proposed by Cornwall (1982) (for references and applications, see j4|). 
This mechanism tames the Landau singularity, but the intrinsic scale m g is an extraneous 
parameter to the approach and has to be determined by other methods, e.g., on the lattice. 
Indeed, the saturated coupling at one loop reads 



the range (500 ± 200) MeV. This scale my be thought of as being the average transverse mo- 
mentum of vacuum partons within the framework of nonlocal quark and gluon condensates. 



with A 2 
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Technically, the role of the IR regulator m g is to "freeze" a s (Q 2 ) at low momenta Q 2 < A 2 , 
so that the coupling ceases to increase and flattens out. 

An alternative approach to protect perturbatively calculated observables was developed 
by Sterman and collaborators (references and applications in Q), in which a s singularities 
were screened by Sudakov-type damping exponentials e~ s that encapsulate the effects of 
gluonic radiative corrections. These are resumed contributions stemming from two-particle 
reducible diagrams (giving rise to leading double logarithms), while two-particle irreducible 
diagrams (entailing subleading single logarithms) are absorbed into the hard-scattering part 
of the process. Because e~ 5 drops to zero faster than any power of In (/i 2 /AQ CD ), this 
provides sufficient IR protection against Landau singularities both for mesonic amplitudes 
as well as for baryonic ones (see [4j for a review). Indeed, in the axial gauge, all Sudakov 
contributions due to unintegrated transverse momenta of the gluon propagators exponentiate 
into suppressing Sudakov factors that suffice to render the Landau singularity of the one- 
loop running coupling in the form factors of the pion and the nucleon harmless, though 
the latter case is theoretically much more involved — see pj]. The main advantage of this 
approach relative to the previous one is that it contains an intrinsic IR regulator in terms 
of the impact parameter characterizing the inter-quark separations in the hadrons. Hence, 
there is no need here for external IR regulators (like m g ) because IR protection is provided 
in situ by the Sudakov effect, the underlying idea being that gluons with wavelengths larger 
than the impact parameter "see" the quark configuration as a whole, i.e., in a color-singlet 
state and hence decouple. It is worth noting that the Sudakov suppression can be simulated 
by a saturated running coupling of the form of Eq. (j7j), if we assume that the IR cutoff 
parameter (i.e., the effective gluon mass) has to be adjusted to each individual value of the 
varying external momentum — for a mathematical version of this idea see Q. Let us close 
this discussion by remarking that the application of this type of approach to the timelike 
region using analytic continuation is not straightforward owing to the fact that not only the 
running coupling has to be analytically continued but also the whole Sudakov form factor. 

A bridge here to what is to come in the next section with regard to analytization. In 
separate parallel developments (see [E| for references), Krasnikov and Pivovarov (1982), and 
independently Radyushkin (1982), have obtained analytic expressions for the one-loop run- 
ning coupling (and its powers) directly in Minkowski space using an integral transformation 
from the spacelike to the timelike region reverse to that for the Adler D-function. In both 
approaches a resummed A-parametrization for a s (Q 2 ) in the spacelike region was employed 
in order to construct in the timelike region an expansion for R(q 2 ) in which all (tt 2 /L 2 ) n - 
terms (L = lns/Aq CD ) are summed explicitly. This way, an analytic coupling in the timelike 
region was derived. This sort of analytic coupling in Minkowski space was rediscovered in the 
context of the resummation of fermion bubbles by Ball, Beneke, and Braun (1995), partly in 
connection with techniques and applications to the r hadronic width (references are given in 
(H). Still other approaches to avoid ghost singularities in the running coupling are referenced 
in 0,8- 



IV. INTO ANALYTIZATION 

Shirkov and Solovtsov (1996) 0] have invented an analytic coupling based solely on 
Renormalization-Group (RG) invariance and causality (spectrality) in terms of a disper- 
sion relation, generalizing the earlier work by Bogoliubov, Logunov, and Shirkov for QED 
(the reference can be found in [6j). The key ingredient of their framework is the spectral 
density p/(cr) with the aid of which an analytic running coupling can be defined in the 
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Euclidean region (— q 2 = Q 2 ) using a Kallen-Lehmann spectral representation: 



The same spectral density p/(cr) = Im[/(— cr)}/n is then used to define the running coupling 
also in the timelike region by means of a dispersion relation for the Adler function. 
At the one-loop level, one obtains in the Euclidean region [6| 



while its counterpart in Minkowski space reads 



StiCs) = / '^-da = - arccos —. Ls (10) 



where L = In (Q 2 /Aq CD ) and L s = In (s/Aq CD ) . Both couplings are analytic functions of 
their arguments. The analyticity of A\(Q 2 ) is ensured by the second power-behaved term 
in (Q which removes the Landau (ghost) pole, whereas 2li(s) does not contain any singular 
term at all. This analytization concept was used earlier in QED by Redmond (1958) (see 
0] for the reference) in connection with the elimination of ghosts in propagators, based on a 
term-by-term perturbation expansion in the coupling constant. The graphical representation 
of the removal of the Landau ghost in the coupling is shown in Fig. [TJ 




FIG. 1: The Landau singularity, marked by the symbol X in the left panel, is absent by construction 
in the analytic approach of Shirkov and Solovtsov (right panel). 



V. THE BIRTH OF APT. FROM A CONCEPT TO A PARADIGM 



The simple analytization concept for the one-loop running coupling in the spacelike region 
has been extended by its inventors and their collaborators to higher loops and has been 
applied to several characteristic hadronic processes. Examples are the inclusive decay of 
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a r-lepton into hadrons, the momentum-scale and scheme dependence of the Bjorken and 
the Gross-Llewellyn Smith sum rule, T decays into hadrons, et cetera (the corresponding 
references can be found in In Fig. [5J we show the first power (left panel) and the 

second power (right panel) of the analytic couplings in the Euclidean (calligraphic notation) 
and in the Minkowski (Gothic notation) space. All couplings are evaluated at the one-loop 
level. One appreciates from this figure that the spacelike and timelike analytic couplings 




FIG. 2: Analytic couplings in APT: spacelike region — A n (Q 2 ) and timelike region — 2l n (s). The 
left panel shows the first power (index), while the second power (index) is given in the right panel. 

have a distorted symmetry for finite values of their arguments, though asymptotically, i.e., 
for L,L S — > oo, they tend to the conventional coupling: A„(L — > oo) =^> 2l n (L s — > oo) =>- 
a n (L — > oo) for n G N, bearing in mind that a = (b /4n) a s . As we will see later, this 
property is valid also for real indices v G M. 

Meanwhile, a systematic approach — a new paradigm for perturbative calculations in 
QCD — has emerged, termed Analytic Perturbation Theory (APT). In a nutshell, it can 
be summarized as follows. The standard and the analytic couplings at one loop can be 
defined jH| recursively for any integer power (index) nGNby means of 





(n-l)\ V dk 
where 

• a n (standard QCD with a = l/L, n G N: power) 

• A n (analytic in spacelike region, n G N: index); k = L = ln(Q 2 /A 2 ) 

• 2t n (analytic in timelike region, n6N: index); k = L s = ln(s/A 2 ) . 
On the other hand, in terms of the spectral density 
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one has 
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where F(z, n) is the transcendental Lerch function serving in (14a) as a Landau-pole remover 
(Bj. In APT, hadronic quantities defined (in the Minkowski region) by § f(z)R(z)dz with 

R PT (Z) = dmOC{z) ^ TZ APT (Z) = dm%m(z) (14) 

n n 

have no Landau singularities in the Euclidean region, because the spacelike couplings A n (z) 
are analytic functions. One appreciates that APT transforms a power-series expansion into 
a non-power-series, i.e., a functional expansion. 

Let us close this section by making some important remarks (for references, see (BJ): 



• Analytization in Euclidean space means subtraction of the Landau pole (at one loop). 

• Analytization in Minkowski space amounts to summation of 7r 2 terms. 

• Observables become non-power-series expansions (i.e., n is an index not a power): 
T>(Q 2 ) = J2 n d n -A n (Q 2 ) (Euclidean space), 7Z(s) = J2 n d n ^- n (s) (Minkowski space). 

• The elimination of ghost singularities in APT appears as a result of causality (spec- 
iality) and RG invariance, i.e., the pole remover is not introduced by hand. 

• Two-loop expressions for the analytic couplings are possible by means of the Lambert 
function (Magradze (2000)). 

• Still higher orders can be obtained with an approximate spectral density and employing 
numerical integration (Shirkov (1999)). 

• Incorporation of heavy-quark flavors possible ("global APT" — Shirkov (2001)). 

• Modifications of APT in the deep IR are possible via a spectral density which contains 
(nonperturbative) power corrections (Alekseev (2006); Nesterenko and Papavassiliou 
(2005); Cvetic and Valenzuela (2005)). 

VI. MORE SCALES, MORE RIDDLES 

The above exposition suffices to demonstrate the usefulness of APT. But where are its 
limitations? Can the new paradigm be applied to more complicated cases that contain more 
than one large momentum scale? And can APT accommodate the crucial effect of evolution? 
As we cannot sidestep these problems, we will have to dig deeper into analytization. This is 
not a terminal difficulty of the analytic approach — and we can even benefit from it because 
doing so, we will ultimately expand APT to noninteger powers of the coupling. It is to these 
questions to which we focus now our attention. 
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To start with, we consider as a "theoretical laboratory" the factorized part of the spacelike 
pion's electromagnetic form factor. In NLO perturbation theory, this quantity can be written 
as a convolution [A ® B = f dzA(z)B(z)] of a hard-scattering amplitude T and two pion 
distribution amplitudes, ip™, <p° nt , representing the incoming and outgoing pion bound states 
of quarks with longitudinal momentum fractions x and y, respectively. T is the amplitude 
for a collinear valence quark-antiquark pair, struck by a highly virtual photon with the large 
momentum Q 2 , and contains hard-gluon exchanges characterized by the strong coupling a s 
to the considered order of expansion. Then, we have 



F? ct (Q 2 ) = <ft(x,,4) ®T£ LO (x.y.g 8 ;^,/^) ® vTfati) > (15) 

where 

(p n (x, fx 2 ) = 6x(l - x) 1 + a 2 (ju 2 ) C 2 3/2 (2x - 1) + a 4 (/U 2 ) Cf /2 (2x - 1) + . . . (16) 

is the leading twist- 2 pion distribution amplitude containing the nonperturbative informa- 
tion on the pion quark structure in terms of the Gegenbauer coefficients a n at the initial 
scale fx 2 ~ 1 GeV 2 , and T is the sum of all quark-gluon subprocesses up to the order a 2 
(its explicit form can be found in the second paper in [2j and also in [3]). Note that T 
depends explicitly on a s , whereas ip n has a more complicated dependence controlled by the 
evolution equation. Besides, beyond the LO of perturbation theory, T depends not only on 
the large external momentum Q 2 , but also on a second auxiliary momentum scale fip which 
separates the hard (perturbative) from the soft (nonperturbative) dynamics on account of 
a factorization theorem. Invariably, the soft parts, i.e., the pion distribution amplitudes, 
have to be evolved with appropriate anomalous dimensions from the initial scale fj 2 to the 
actual scale of observation Q 2 . As we will discuss in the next section, this additional scale 
has influence on the analytization procedure that cannot be covered by APT. 

The first application of APT to Fj act (Q 2 ) was considered in [21, in which the strong; 
coupling and its powers were replaced by their corresponding analytic images: a n (L) — > 
(Ai(L)) n (Naive analytization [3j). This treatment is strictly speaking incorrect because 
[*4.i(L)] n ^ [a™(L)]An, but phenomenologically it works rather good An improvement 
of this approach was presented in where a Maximal analytization of this quantity was 
performed. This procedure associates to the powers of the running coupling their own dis- 
persive images, i.e., [a™(L)]M ax -An = A n (L). It is worth emphasizing that the analytization 
treatment was applied not only to the hard part T but also to the anomalous dimensions, 
governing the evolution of the pion distribution amplitudes, which can be computed order- 
by-order in perturbation theory. The differences between the two analytization procedures 
becomes apparent by comparing the following two expressions: 

Naive Analytization 

[Q 2 T n (x, y, Q 2 ; /4 , ArQ 2 )] Naivc __ An = A? (A R Q 2 ) 4° } (*, y) 

~A?\\ R Q 2 ) 



+ 



^Uv^kM) (17) 



4vr " a V^'^'Q 2 

Maximal Analytization 

[Q 2 T H (*, y, Q 2 ; /4 , ArQ 2 )] Max An = A? (X R Q 2 ) $ (x, y) 

Af{KQ 2 ) (l) ( » 2 \ 



9 




0.1- ; . 1 



10 20 30 40 50 10 20 30 40 50 

Q 2 [GeV 2 ] Q 2 [GeV 2 ] 

FIG. 3: Results for Q 2 Fj act (Q 2 ) vs. Q 2 with /4 = Q 2 , $ = 5.76 GeV 2 and using Naive Analyti- 
zation (left panel) and with Maximal Analytization (right panel) for various renormalization-scale 
settings described in Q]. All predictions have been calculated with the Bakulev-Mikhailov-Stefanis 
pion distribution amplitude (see 0] for references). 



where ArQ 2 = /j r serves as a renormalization scale and /x F denotes the factorization scale of 
the process. The explicit form of the amplitudes and is irrelevant for our arguments 
here; it can be found in Furthermore, in [2j , the (naive) analytization procedure was 

also applied after exponentiation to the Sudakov form factor in next-to-leading logarithmic 
approximation of the cusp anomalous dimension. 

Some technical comments: (a) The pole remover in the one-loop running coupling does 
not change the leading double logarithmic behavior of the Sudakov factor, though the cusp 
anomalous dimension in the Sudakov exponent receives additional subleading contributions, 
(b) Due to the absence of ghost singularities, the Sudakov exponents are analytic functions 
of Q 2 . (c) The perturbatively calculable part of the pion's electromagnetic form factor 
Q 2 F^ act (Q 2 ) receives the bulk of its contributions from small transverse distances b below 
b ■ Aq CD < 0.5, i.e., it saturates in exactly that region, where the application of QCD 
perturbation theory is valid 0. However, the full treatment of the Sudakov-type gluon 
resummation in impact space faces special subtleties which have not been fully resolved 
until now. First results within a toy model in momentum space (see Appendix C of 0) 
seem to indicate a conflict between analytization and exponentiation of soft gluons. 

Phenomenologically, already the "naive" analytization procedure and even more the 
"maximal" application of APT have important consequences: (i) the artificial increase of the 
spacelike electromagnetic pion form factor towards low Q 2 values — caused by the Landau 
singularity — is absent, (ii) the sensitivity of the results on the choice of the renormalization 
scheme and the scale setting adopted is extremely reduced relative to the standard power- 
series expansion of perturbative QCD [Shirkov and Solovtsov (1998) and This behavior is 
shown graphically in Fig. [3] for the scaled expression of Fj act (Q 2 ) employing these two analy- 
tization procedures and including leading-order evolution of the pion distribution amplitude 
in APT. The predictions shown were calculated for different renormalization-scale settings, 
including the "default" choice, /i R = Q 2 (dashed line), the Brodsky-Lepage-Mackenzie scale 
setting (dotted line), the so-called ay scheme (dash-dotted line), et cetera, and employing the 
Bakulev-Mikhailov-Stefanis pion distribution amplitude which has ^2(^0 = 1 GeV 2 ) = 0.20, 
a 4 (/iQ = 1 GeV 2 ) = —0.14, a n (nl = 1 GeV 2 ) ^ (n > 2). All these issues are discussed 
in detail in [3|. The main observation from this figure is that applying APT in its maximal 
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form to a typical exclusive quantity, like F^ act (Q 2 ), not only provides IR stability, it also di- 
minishes its renormalization scheme and scale-setting dependence already at the NLO level. 
Indeed, the predictions obtained in different schemes do not differ significantly from each 
other in a wide range of momenta from low to large Q 2 values. 

To further appreciate the advantages of the analytic scheme relative to the standard 
power-series expansion, bear in mind our statements in Sec. [TT] in the language of infor- 
mation theory. Within this context, one may consider the variation of the predictions — the 
uncertainty spread — obtained with different renormalization schemes, as "noise", the predic- 
tion for Fj act (Q 2 ) itself being the "signal" . It becomes obvious from the right panel of Fig. [3] 
that maximal analytization within APT provides a much better dynamical range of quality 
because the signal-to-noise ratio is strongly enhanced in the whole Q 2 range, especially at 
low Q 2 . 



VII. FRACTION ALIZING APT 

In the previous section we have discussed the advantages of APT relative to the standard 
QCD perturbation theory and considered a typical exclusive quantity, notably, the factoriz- 
able part of the pion's electromagnetic form factor that exhibits a dependence on potentially 
two large momentum scales: Q 2 and the factorization scale /xp (that may also serve as the 
evolution parameter). Until now our discussion has not addressed the factorization-scale de- 
pendence under the proviso of analytization. But the considered example makes it obvious 
that, ultimately, one has to accommodate the analytization of terms like 



Z[L] = e rlL] 7$> da -> [a s (L)] 7o/2A) < — > RG at one loop 

[a s (L)] n In [a s (L)] < — > RG at two loops 

[a s (L)] n L m < — > Factorization 

exp [— a s (L)F(x)] < — > Sudakov resummation 



(19) 



typically appearing in perturbative calculations beyond LO, as in our example. Though 
such terms do not modify ghost singularities, say, the position of the Landau pole in the 
one-loop approximation, they do contribute to the spectral density. To include such terms 
into the dispersion integral, one is actually forced to modify the analyticity requirement and 
demand not the analyticity of the running coupling and its powers but rather the analyticity 
of the quark-gluon amplitude as a whole [{| . This generalized encompassing version of the 
analyticity requirement demands that all terms that may contribute to the spectral density, 
i.e., affect the discontinuity across the cut along the negative real axis — oo < Q 2 < 0, must 
be included into the analytization procedure 0, H, 0] . The use of this principle allows the 
inclusion into the dispersion relation of logarithmic terms, of the sort ln(Q 2 / /j p ), or products 
of such logarithms with powers of the running coupling. It turns out that precisely such terms 
are tantamount to fractional (in fact, real) powers of the strong coupling, pertaining to FAPT 
@,|. 

Before studying the consequences of this generalized analyticity requirement for observ- 
ables, let us first present the different analytization concepts, we have discussed above, in 
mutual comparison — Fig. HI In this figure, the linear operations Ap and Am define, respec- 
tively, the analytic running couplings in the Euclidean (spacelike) 

r°° Jf)(fj\ 

AeK)]=^ ) with A n(Q 2 ) = / ^^Mda (20) 

Jo °~ + V 
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FIG. 4: Implementation of analyticity in standard QCD perturbation theory (PT), APT (a), and 
in FAPT (b) and (c). 

and the Minkowski (timelike) region 

Am [afo] = 2t« with Zg>(s) = [°° ^ da , (21) 

J s °" 

where the loop order of the expansion is denoted by the superscript (/). The above analyti- 
zation operations can be represented by the following two integral transformations from the 
timelike region to the spacelike one: 

D[X®]=JW with ^^f* (22) 

Jo [a + Q 2 ) 

and for the inverse transformation 

R[J®]=V® with aW (a)= 1 r +l£ ^l d a. (23) 

2m J_ s _ i£ a 

These two integral transformations are interrelated: 

DR = RD = 1 (24) 

for the whole set of analytic images of the powers of the coupling in the Euclidean, as well 
as in the Minkowski space, {^4 n ,2l n }> respectively, and at any desired loop order of the 
perturbative expansion. Then, Eq. (14a) generalizes from integer indices n to any real index 
i/6lto read (Euclidean space) 

A (L) - 1 - F{6 ~ L ' 1 ~ v) (25) 

where the first term corresponds to the conventional term of perturbative QCD (with L = 
ln(Q 2 /A 2 )) and the second one is entailed by the pole remover (cf. l/(e L — 1) at one loop). 
The function A U (L) is an entire function in the index v and has the following properties: 

• ML) = 1- 

• A-m(L) = L m = a~ m (L) for m G N. Hence the inverse powers of the analytic coupling 
and the original running coupling coincide). 
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• A«(±oo) = for m > 2, m G N. 

. A V (L) = - [1/I»] Er=o C(l - f - r) [(-L)7r!] for |L| < 2vr. 

Note that the reduced transcendental Lerch function F(z, v) is related to the Lerch function 
i/, m) via z$(z, z/, 1) = F(z, v). For a positive integer index, v = m > 2, one has the 
relation 



F(z,l 

so that substituting Eq. (J25|) in 



I -.1 



one arrives at 

Ar.iL) = {-IT A 



-L) . 



(26) 



(27) 



In Minkowski space, the analytic images of the running coupling are completely deter- 
mined by elementary functions |8j (L s = ln(s/A 2 )): 



21, (L a ) 



sm 



1) arccos ( L 



7H + L2 



tt(z/ - 1) (vr 2 + L 2 ) 



(1/-1J/2 



(28) 



from which, for example, we get 2lo(L s ) = 1, 2l_i(L s ) = L s ; 2l_2(-^) = L 2 — 2l m (L) = 
(-l) m 2l m (-L) for m > 2; 2l m (-oo) = An(-oo) = 5 m ,i; 2t m (oo) = An(oo) = for m G N. 
For the value L = L s = 0, the analytic couplings in both regions are interrelated by the 
equation 



A(o) 



i)CM 



21,(0) , 



where 



21,(0) 



sin [(i/- 1)tt/2] 



1) 7T y 



(29) 



(30) 



and where ((u) is the Riemann ( function, with the coefficient in the bracket providing a 
quantitative measure for the magnitude of the distortion of the "mirror symmetry", men- 
tioned before, for any v G R. A graphical illustration of the analytic couplings A v and 21, is 
given in Fig. [5] which shows the rate of change of these functions with respect to the index 
v and the argument L (or L s ). 

The major characteristics of FAPT in comparison with the standard perturbative expan- 
sion in QCD and the original APT are collected in Table HJ while further details can be found 

m Baa. 

Up to now we have actually concentrated our attention on the one-loop approximation of 
the running coupling and its analytic images. But the index derivative in FAPT gives us the 
possibility to obtain two-loop expressions from those at one loop. To achieve this goal, let us 
recall the two- loop RG equation for the standard QCD normalized coupling am = b a^ /Arc 



da(2)(L) 
dL 



~ a (2)(£) [l+CiCJ(2)(L)] 



which assumes the following functional form 

1 



1(2) 



+ Ci In 



1(2) 



1 + CiO( 2 ) 



with c\ = 



L . 



(31) 



(32) 
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FIG. 5: Comparison of the Euclidean (left panel) and the Minkowski (right panel) analytic cou- 
plings, A v {k = L) and %l u (k = L s ), respectively, for incremental changes of the index v in the range 
2 to 3. 



TABLE I: FAPT versus standard QCD perturbation theory (SPT) and APT. 



Theory 


SPT 


APT 


FAPT 


Space 




{A m , 2lm} meN 


{A u , 2ly}j, gR 


Series expansion 


Ef m a m (L) 

m 


EfmA m (L) 

m 


fm A m (L) 

m 


Inverse powers 


[a(L)]- m 




A-m(L) = L m 


Multiplication 


a»a v = a^ +v 






Index derivative 


a v ln fc a 







This equation has the following exact solution (Magradze 2000) 

tt(2)(L) = i + N»(t)) ' (33) 

where zw(L) = (1/ci) exp (— 1 + in — L/ci) with Wk (k = 0, ±1,...) being the Lambert 
function, defined by z = W(z) exp (W(z)), where its different branches are depicted by the 
index k. 

With the above expression in our hands, we can immediately obtain [8[ the spectral 

/r)\ 

density at the two- loop level, pi (cr), by virtue of Eq. (}T2"]) . viz., 

p^(a) = -Im[a^(L-zn)] . (34) 
Using this expression in Eqs. fl20l and fl2Tj) . we have 




from which we derive {T> = d k /dv k ) 

A? (L) = A? + Cl V A u % + cj{V* + V- 1) A[l, + O (X> 3 Al%) , (36) 
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where we have expanded a(o\{L) in terms of am(L). The above expression can be generalized 
to any real index v, in both the Euclidean and the Minkowski space, to read 



A {2) 
2l< 2) 



A {1) 



21 



(i) 



c\vT> 




u + l 
2! 



V + V — 1 



A {1) 
2t (1) 



+ O (c?) 



(37) 

Though these expressions are not exact, the achieved accuracy is of the order of about 
1% down to the value L = L s = and, hence, they are sufficient for most practical QCD 
applications. A more accurate expression for the two-loop Minkowski coupling which includes 

the contribution of the next higher order O (^D 4 21^4 J is given in [5]. Moreover, one can 

find there the analogous three-loop expression with a similarly good quality of convergence. 

The particular advantage of FAPT is the possibility to perform the analytization 
of powers of the running coupling multiplied by logarithms of L and its powers, e.g., 
[(ct(2)(L)Y £ m ] an = £u)n(L). As we emphasized before, such terms appear in NLO cal- 
culations of typical hadronic amplitudes, like form factors, or when taking into account 
evolution effects. Table [Til displays in approximate form the analytization of such expres- 
sions. Their usefulness will become evident below, when we will consider the influence of 
such terms in the calculation of the factorized part of the pion's electromagnetic form factor. 



VIII. FAPT APPLICATIONS 



1. Having set up the calculational tools of FAPT, we now detail its concrete application, 
starting with an example in the spacelike region. Consider again Eq. ([18]) — this time includ- 
ing into the analytization procedure also the logarithmic term (labeled by the acronym KS 
0): ln(Q7/4) = ln(A R Q 2 /A 2 ) - ln(A R /4/A 2 ). Then, we find (x = 1 - x) 



[Q 2 T n (x,y,Q 2 ;^X R Q' 



2M An 



+ A^ArQ 2 ) fc Ft g) (a . y) ( 6 + 21n(^)) 



12W(0), 



4 2) (ArQ 2 ) + (d 



47T 



/4 



47T 



with 



A?> (Q 2 ) = 4 2) (Q 2 ) - 4 2) (Q 2 ) In [Q 2 /A 2 ] 



,(2) 



(2) 



(3J 



(39) 



encoding the deviation from the result obtained with the Maximal Analytization procedure 
and where 



Cf (Q 2 )i 



K 2) (Q 2 )n»(f 



An 



47T 



Jks u o 



a {2) (Q 2 



(Q 2 ) 



An 



KS 



(40) 



Carrying out the KS analytization, we get the approximate expression 



£^ pr(2) (Q 2 ) 



47T 



47T 



4 2j (Q 2 ) +ciirfc (Q 2 ) 

On 



(41) 
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TABLE II: Calculational rules for FAPT with L = In (Q 2 /A 2 ), m G N, and v G R. J x is a Bessel 
function of the first kind and ip(z) = dlnT(z)/dz denotes the Digamma function. 



Standard QCD PT 



QCD FAPT 



L" 



fy(L) ln m [a (l) (L)] 



a v {2) {L)L 



a<( 2) (L)L 2 



exp [-xa(L)] 



A { u ] {L)=Al L '(L) + c 1 uVA K ;i 1 (L)+c{ 



(i) 



3 /I CO 



(*±l)«r>2 



£> 2 + I? - 1 



+ [V 6 A 



v+3 



45 (L) = ^(L) + Cl D^(I) + (P 2 ^ 



l (2) 



(2)/ 



^2 ,(2) 



= ^ 2 (L) + 2 d IM^L) + c 2 V 2 A K u'(L) - 2c{ A K u'(L) 



I (2) 



I (2) 



2-n2 ,(2) ( 



,2 ,(2), 



^ 2 (L)+p l i/0^£ 1 (£)+^ 



+ ViE m =ie- mi ^^ forL>0 



(which can be further improved by substituting the ci-term in Eq. (HIT) by its two-loop 
expression : C™l' m [L] = Af ] [L] + a V Af [L] —see Q), where [cf. Table [HI a^ 2 (L)L] 

n>0 

and is the Riemann Zeta-function. The main effect of including the logarithmic term 
into the analytization procedure of F^ act (Q 2 ) is to suppress its sensitivity on the variation 
of the factorization scale to a minimum. It was shown in [7] that varying the factorization 
scale from 1 GeV 2 to 10 GeV 2 , the form factor changes by a mere 1.5 percent reaching 
for a (hypothetical) factorization scale of 50 GeV 2 just the level of about 2.5 percent. To 
conclude, using FAPT, the dependence of a perturbatively calculated QCD quantity on all 
perturbative scheme and scale settings is diminished already at the NLO level. 

2. The second application deals with the Higgs boson decay into a bb pair in the timelike 
region. We will present here only the mains results and refer for the technical details to (EJ . 
A dedicated analysis of the uncertainties involved in the calculation of the decay width of 
the Higgs boson into bottom quarks was recently given by Kataev and Kim — second entry 
ofRef. Q. 



V(2)C(-n-l)- 



dn 



[-m(Q 2 /A 2 )]' 

r(n + i) 



(42) 
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Consider now the decay of a scalar Higgs boson to a bb pair in terms of the quantity 
Rs at the four- loop level from which one can then obtain the width r(H — > bb). Note 
that the application of FAPT in Minkowski space is not plagued by Landau singularities. 
However, the analytic continuation from the spacelike to the timelike region induces so- 
called "kinematical" ir 2 terms that may be comparable in magnitude with the higher-order 
expansion coefficients and have, therefore, to be summed. The first step in our approach is 
the correlator of two scalar currents Jj? = ^b^b for bottom quarks with mass rnj,, coupled 
to the scalar Higgs boson with mass Mr, and where Q 2 = —q 2 , i.e., 



H(Q 2 ) = (4vr) 2 * / dxe^ x {0\ T[ J!(x)J^(0)] |0> 



(43) 



Then, R s (s) = ImU( 



-s — ie)/(2ir s) and one can express the decay width in terms of R$: 



r(H -> bb) 



Ay/2 



-M n ml(M n )R s (s 



Ml) 



(44) 



7T 



The quantity Ro, is obtained from the Adler function D via analytic continuation into the 
Euclidean space using the transformation Am (or, equivalently, the integral transformation 
R). The reason for this detour is that in the Euclidean region perturbation theory works. 
Hence, we write 



D S (Q 2 ^ 2 ) = 3m 2 (Q 2 



l + $> n (Q 2 //i 2 ) <(/i 2 ^ 



n>l 



(45) 



The second step is to calculate 



R s (s) = R S (Q 2 = s^ 2 = s) = 3m 2 b (s) [l + (46) 

n>l 

using the techniques developed by Gorishnii et al. (1990) and Chetyrkin et al. (1997) fu\ . 

As already mentioned, the coefficients r n contain characteristic tt 2 terms originating from 
the integral transformation R of the powers of the logarithms entering D$. The latter is 
related to Rs(s, s) by means of a dispersion relation. Notice that these logarithms have two 
different sources: one is the running of a s in D§, while the other is related to the evolution 
of the heavy-quark mass m 2 (Q 2 ). As a result, the coefficients r n in (j4"6"|) are connected 
to the coefficients d n in Da, (calculable in Euclidean space) and to a combination of the 
mass anomalous dimension 7$ and the /3-function coefficients bj multiplied by ir 2 powers. 
Appealing to the explicit results given in we remark that the total amount of these terms 
is of the order of the coefficient d±. FAPT can by construction account for these terms to all 
orders of the perturbative expansion. This is because in FAPT one does not need to expand 
the renormalization factors into a truncated series of logarithms; instead one can transform 
them by means of the AM-operation " whole" . 

The running mass in the /-loop approximation, m(j), can be cast in terms of the 
renormalizat ion-group invariant quantity my) to read 

™UQ 2 ) = rh% [a s (Q 2 )Y° f {l) (a s (Q 2 )) , (47) 
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where uq = 270/60 and the expansion of f(i)(x) at the three-loop order is given by 



l + a s 



h ( 71 



26 V 61 



To 
60 



+ a 



b\ 



°16b% 



6 61 



To 
60 



Ti 
61 



6062 


^72 


ToV 


&1 


u 





0(a2) • 



(48) 



Expanding the running mass in a power series according to 



in 



(0 



S)W ) 

and setting /i 2 = Q 2 , we find 

[3m 2 ]-;i)«(g 2 ) = (ai')(g 2 )) 

with 



(« S (Q 2 )) ! 



m>l 



(49) 



n>l 



4 WW ° + E A £^W)) 



m>l 



min[/,m— 1] 

e (0 + V tLe (Z) 
fc>l 



(50) 



(51) 



Note that we have separated the mass-evolution multiloop effects (third term of Eq. (1501) ) 
from the original series expansion of D (truncated at n = I). These contributions are 
represented by the second term on the RHS of Eq. (I501) . In practice, for Q > 2 GeV, i.e., 
for a s < 0.4, the truncation at m = I + 4 of the summation in fj4*9l is sufficient, given that 
the truncation error is less than 1 percent. 

By applying the analytization operation A M to the quantity Ds\Q 2 ), we finally obtain 



(i)MFAPT 



A M [D 



3m 2 (l) 



(0i 
s J 



n>l 



m>l 



(52) 



where we have used the shorthand notation 



cms 



o!!>W 



5;' apW ' 



(0 / Jih 



(53) 



and the couplings C^ +UQ , all 



This expression contains, by means of the coefficients A 
renormalization-group terms contributing to this order. On the other hand, the resummed 
7r 2 terms are integral parts of the analytic couplings a^ +vo . 

The results for R S (M^) as a function of the Higgs mass Mh, calculated with different 
perturbative approaches in the MS scheme, are shown in Fig . El The long-dashed curve 
shows the predictions obtained in the second paper of Ref. HI] using standard perturbative 
QCD at the I = 4 loop level of the expansion. The solid curve next to it represents the 
FAPT prediction (cf. fl52|) ). including in the second sum all evolution effects up to m = I + 4 
and fixing the active flavor number to Nf = 5. On the other hand, the contributions of 
the higher-loop renormalization-group dependent terms are accumulated in the coefficients 
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i?s(M H ) [GeV : 



32 



30 



2 8 ■ 



26 




32 



30 




60 80 100 120 



60 80 100 120 



FIG. 6: Predictions for the quantity Rs(M^) calculated in different perturbative approaches 
within the MS scheme: Standard perturbative QCD 11] at the loop level I = 4 (dashed red line 



with An /= 5 = 231 MeV), Broadhurst-Kataev-Maxwell ("naive non-Abelianization" ) including the 
O({a s ) vo A 4 (a s ))-terms, [1]— (dotted green line with Ajv /=5 = 111 MeV), and FAPT [5] obtained 
from Eq. ([52]) for Nf = 5 (solid blue line), displayed for two different loop orders: I = 2 (left panel, 
A 



Nf 



263 MeV) and I = 3 (right panel with A 



N r 



261 MeV). The value of Ajv f=5 MeV in 



all cases corresponds to 21^ {s = rn 2 z ; Nf = 5) = 0.120. 



2 . 



^1 by means of the parameters ji and by Obviously, the standard perturbative QCD 
approach and FAPT yield similar predictions for this observable, starting with the two- 
loop running. The reason for the slightly larger FAPT prediction lies in the fact that the 
coefficients a u contain by means of the index uq the resummed contribution of an infinite 
series of 7r 2 -terms (as well as all 70 and &o terms) that renders them ultimately smaller than 
the corresponding powers of the standard coupling. Finally, the dotted green curve about 
8% below the other lines, in both panels of Fig. [HI illustrates the estimate obtained with the 
"naive non-Abelianization" and an optimized power-series expansion that makes use of the 
"contour integration technique" [lj]. 



IX. CONCLUSIONS 



FAPT derives from APT and generalizes it both conceptually and technically. Conceptu- 
ally, it extends the analyticity requirement from the running coupling and its powers to the 
hadronic amplitude as a whole. Technically, it extends the formalism from integer powers 
of the coupling to real ones in the spacelike as well as in the timelike region. More recently, 
progress has been achieved by endowing FAPT with a varying flavor number across heavy- 
quark thresholds: "flavor-corrected global FAPT" (see [To| for explanations and details). Yet 
FAPT is not a complete formalism. We still do not really know how to consider Sudakov 
gluon resummation within FAPT because under the analytization requirement resummation 
does not mean exponentiation. This important issue has yet to be understood and devel- 
oped. Also the application of FAPT to compute power corrections to the spacelike pion's 
LO electromagnetic form factor and the inclusive Drell-Yan cross section, started in [9J , has 
so far not been exploited beyond the leading power and for other hadronic quantities. 
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